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ABSTRACT 

Transiting extrasolar planets constitute only a small fraction of the range of stellar 

■ systems found to display periodic, shallow dimniings in wide-field surveys employ- 
rS \ ing small-aperture camera arrays. Here we present an efficient selection strategy for 

■ follow-up observations, derived from analysis of the light curves of a sample of 67 Su- 
per WASP targets that passed the selection tests we used in earlier papers, but which 
have subsequently been identified either as planet hosts or as astrophysical false posi- 
tives. We determine the system parameters using Markov-chain Monte Carlo analysis 
of the SuperWASP light curves. We use a constrained optimisation of combined 
with a Bayesian prior based on the main-sequence mass and radius expected from the 
2MASS J — H colour. The Bayesian nature of the analysis allows us to quantify both 
the departure of the host star from the main-sequence mass-radius relation and the 
probability that the companion radius is less than 1.5 Jupiter radii. When augmented 
by direct light curve analyses that detect binaries with unequal primary and secondary 
eclipses, and objects with aperture blends that are resolved by SuperWASP, we find 
that only 13 of the original 67 stars, including the three known planets in the sample, 
would qualify for follow-up. This suggests that planet discovery "hit rates" better than 
one-in-five should be achievable. In addition, the stellar binaries that qualify are likely 
to have astrophysically interesting stellar or sub-stellar secondaries. 

Key words: methods: data analysis - stars: planetary systems - techniques: photo- 
metric 
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1 INTRODUCTION 

Wide-field photometric surveys using commercial, small- 
aperture camera lenses and large-format CCDs have 
yielded seve ral discoveries of new transiting extraso- 



lar planets (lAlonso et all l2004l: iMcCuUough et all 
O'Donovan et alj bOOeT 1 Collier Cameron et al. 



Bakos et al 



2007al: iBakos et alJ l2007bl : iBurke et all 



200e 



2007 



2007 



O'Donovan et al.l 



20071 1 ~ Deeper surveys based on light 
curves from the Optical Gravitational Lensing Experiment 

i OGLE) projec t have produced a further fi v e such planets 
Konacki et all l2003l: iBouchv et alJ l2004l . iKonacki et alJ 
l2004l : |Pont et all 12004 iKonacki et al1l2005l) . 

The most efficient foUowup strategies (in terms of tele- 
scope time expended per false positive) are high-resolution 
radial velocity su rveys using 2-m and larger telescopes. 
IPont et al.l l|2005al ) points out that a single spectrum quickly 
unmasks double-lined spectroscopic binaries and rapid ro- 
tation resulting from tidal locking of the primary's rota- 
tion by a stellar-mass companion. It also reveals pressure- 
broadening of strong absorption- line wings, distinguishing K 
giants from K dwarfs. Second-epoch spectra then eliminate 
the narrow, single-lined stellar spectroscopic binaries among 
the surviving main-sequence stars. All the teams engaged in 
such searches repo rt high rates of as trophysical false posi- 
tives. For example. IPont et al.l ()2005a l'l carried out a detailed 
analysis of targets followed up in the OGLE Carina field. 
They found that nearby blended eclipsing binaries, grazing 
equal-mass stellar binaries and transits by planet-sized stars 
near the bottom of the main sequence are the most common 
impostors. 

Several authors have described methods for pre- 
selecting targets for such programmes in order to minimise 
the number of astrophysical fals e positives in samples se- 
lected for radial- velocity foUowup. ISeager fc Mallen-Ornela^ 
l|2003l ) pointed out that the relative radii of the planet and 
star can be estimated from the transit depth, while the dura- 
tions of ingress and egress can break the degeneracy between 
the impact parameter and the star's radius. The transit du- 
ration then yields the radius of the primary, and Kepler's 
third law establishes its density, allowing dwarf-giant sepa- 
ration. The masses of bona-fide main-sequence primaries are 
then estimated from the mass-radius relation. 

In practice, however, the signal-to-noise ratios of tran- 
sit profiles secured with wide-field cameras are seldom good 
enough to determine the duration of ingress and e gress re- 
liably This problem led iTinglev fc SackettI (j2005l ') to pro- 
pose a simpler test, based on transit duration, depth and 
orbital period, for the selection of transit candidates wor- 
thy of foUowup. IPrakd ( 2003') and .Sirko fc Paczyhski (2003) 
used ellipsoidal variations in the out-of-transit light curve 
to eliminate stellar binaries. Even with careful pre-selection 
on transit depth, duration, and primary-star colour and el- 
lipsoidal variations using a combination of all these meth- 
ods, the success rate of such foUowup programmes currently 
runs at about one planet per ten or more stars surveyed 
(jPont et al.ll2005al ). To make the most efficient use of fol- 
lowup facilities, more efficient candidate-selection methods 
are clearly needed. Our goal is to perform efficient selection 
using the discovery observations in conjunction with data 
mining of publicly-available databases of photometry and 



proper motions, without having to perform further time- 
consuming observations. 

Here we present a Bayesian approach to identifying tar- 
gets whose light curve morphologies and infrared colours are 
consistent with expectations for main-sequence stars with 
transiting Jupiter-sized companions. In Sections [2] and [3] we 
describe the SuperWASP photometry and foUowup spec- 
troscopy. The spectra, augmented by photometric tests for 
eliminating blends with nearby eclipsing binaries and graz- 
ing stellar binaries with observable secondary eclipses,yield 
spectroscopic and photometric classifications of 67 high- 
priority transit candidates selected from the 2004 Super- 
WASP transit survey data. 

In Section|4]we analyse the light curves of all 67 stars to 
derive posterior probability distributions for their orbital pa- 
rameters and radii via Markov-chain Monte-Carlo analysis. 
In Section [5] we show that the stellar and planetary param- 
eters obtained for WASP-1 and WASP-2 with this method 
compare well with results obtained from high-precision tran- 
sit photometry with larger instruments. In Section [6] we 
develop candidate selection criteria based on the posterior 
probability that the companion has a radius appropriate to a 
planet, and apply them to the SuperWASP sample. An ancil- 
lary method of dwarf-giant separation based on the reduced 
proper motion and the J — H colour is given in Section [T] 



2 OBSERVATIONS 

Among the 109 SuperWASP transit candidates selected for 
potential spectroscopic foUowup, 67 stars were examined in 
sufficient detail to establish their nature. These 67 objects 
are listed in Table [T] together with their spectroscopic and 
photometric classifications. They we re initially selected for 
foUowup usin g criteria d e scribe d by Collier Camero n et al.l 
(.20061. Christiaii et al.l (|2006| ). iStroct ct al, (200^) and 
iLister et al. I l|2007l ). Among them are the thr ee known planet 



host stars XO-1 jMcCuUough etHI 120061) . WASP-1 and 
WASP-2 (jColUer Cameron et all 120071 ). The remaining 64 
were eliminated as planet candidates using the spectroscopic 
and photometric criteria described here and in Section |3l 



2.1 SuperWASP photometry 

The light curves analysed in this paper were secured be- 
tween 2004 May and September using the SuperWASP wide- 
field camera array. The Super WASP instrument an d data 
pipeline are described in full bv IPoUacco et all l|2006l ). Dur- 
ing 2004, the instrument comprised five f/1.8 Canon lenses 
of 200-mm focal length backed by Andor CCD arrays of 
2048 X 2048 13 /xm pixels. Each camera imaged an area of 
sky 7.8° square. The automated observing routine followed 
a raster pattern sweeping from 3.5 hours east to 3.5 hours 
west of the meridian, returning to a given field every 6 min- 
utes on average. The light curves of most targets comprise 
some 3000 observations secured over a period of 100 to 150 
days. 

The ensemble photometry from each camera on each 
night was calibrated and corrected for extinction, colour 
and zero-point terms using networks of local secondary stan- 
dards. Additional patterns of systematic error affecting all 
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Table 1. Classification criteria for 67 SuperWASP transit candidates. The codes in the Tel/Inst column are E = OHP 1.93-m with 
ELODIE, S=OHP 1.93-m with SOPHIE, I = Isaac Newton Telescope with Intermediate-Dispersion Spectrograph, T=Tautenburg 2- 
m Telescope, W=SuperWASP, X = XO project. Spectroscopic classifications denote whether the spectrum is rotationally broadened 
(Gaussian c > 8 km s~^), single or double, upper limits on radial-velocity variability, and BisVar= Line-bisector variability. The 
presence of interstellar sodium or Class III pressure broadening in the wings of strong lines indicates a giant. Photometric classifications 
are AB=Aperture blend, SE2= Secondary eclipse seen when folded on twice best-fit period, EV2=Ellipsoidal variation seen when folded 
on twice best-fit period. The values of the main-sequence (MS) prior and the probability that the companion radius is less than 1.5-Rjup 
yield the selection criteria in the penultimate column; the final selection including photometric rejection criteria is given in the final 
column. The acronyms RPM and MCMC denote reduced proper motion and Markov-chain Monte Carlo respectively. 



ISWASP Tel/Inst Spec. Phot. RPM MS Prior Probability Select? Select? 

Class Class Class 5 = -2 In (Rp < 1.5Rjup) (MCMC (MCMC 

ViA'h,R*) only) & phot) 



Planets: 

J002040.07-f 315923. 7 S Planet (WASP-1) — 

J160211. 83-1-281010.4 X Planet (XO-1) — 

J203054.12-I-062546.4 S Planet (WASP-2) — 



Dwarf 0.98 0.989 TRUE TRUE 

Dwarf 0.04 0.958 TRUE TRUE 

Dwarf 0.84 1.000 TRUE TRUE 



Sp. Binaries: 






— 


J130322.00-f 350525. 4 


I 


SBl 


— 


J133339.36-I-494321.6 


E 


SBl 





J134815.27-I-464811.0 


E 


Broad 


— 


J174118.30-I-383656.3 


S 


SBl 


— 


J175401. 58-1-322112. 6 


E 


SBl 


— 


J181858.42-I-103550.1 


E 


Broad 


— 


J184119.02-I-403008.4 


T 


SBl 


— 


J203704.92-I-191525.1 


T 


SBl 




J210009.75-f 193107.1 


E 


Broad 




J211608.42-fl63220.3 


S 


SB2 




J212855.03-I-075753.5 


I 


SBl 




J223320.44-I-370139.1 


s 


SBl 




J234318.41-I-295556.5 


s 


SBl 




J010151. 11-1-314254.7 


S,T 


SBl broad 




J023445.65-I-251244.0 


T 


SBl 




J031103.19-I-211141.4 


S 


Broad 




J051221.34-f 300634.9 


S 


Multiple 




J115718.66-f 261906.1 


I 


SBl 




J141558.7H-400026.7 


E 


Broad 




J152131.0H-213521.3 


I 


SBl 




J153135.5H-305957.1 


E 


Broad 




J165949.13-I-265346.1 


S 


Broad 




J172549.13-I-502206.4 


E 


Broad 




J173650.17-I-105557.9 


E 


Broad 




J173748.98-I-471348.7 


I 


SBl 




J174327.81-I-582512.7 


E 


Broad 




J175511. 09-1-134731. 5 


E 


Broad 




J175620.84-I-253625.7 


T 


SB2 




J180010.55-f 214510.2 


T 


SBl 




J182620.36-f475902.8 


E 


SB2 




J202824.02-fl92310.2 


E 


Broad 




J203906.39-I-171345.9 


E 


SBl 




J215802.14-I-253006.1 


S 


SB2 




J231533.56-I-232637.5 


T 


SBl 





Dwarf 


3.87 


0.445 


TRUE 


TRUE 


Dwarf 


0.24 


0.992 


TRUE 


TRUE 


Dwarf 


0.01 


0.864 


TRUE 


TRUE 


Dwarf 


0.86 


0.983 


TRUE 


TRUE 


Dwarf 


0.00 


0.995 


TRUE 


TRUE 


Dwarf 


0.26 


0.905 


TRUE 


TRUE 


Dwarf 


0.86 


0.970 


TRUE 


TRUE 


Dwarf 


0.08 


0.999 


TRUE 


TRUE 


Dwarf 


0.30 


0.998 


TRUE 


TRUE 


Dwarf 


1.42 


0.396 


TRUE 


TRUE 


Dwarf 


0.11 


0.312 


TRUE 


TRUE 


Dwarf 


0.22 


0.998 


TRUE 


TRUE 


Dwarf 


0.10 


0.988 


TRUE 


TRUE 


Dwarf 


17.56 


0.970 


FALSE 


FALSE 


Dwarf 


7.37 


0.034 


FALSE 


FALSE 


Dwarf 


10.52 


0.000 


FALSE 


FALSE 


Dwarf 


7.36 


0.001 


FALSE 


FALSE 


Dwarf 


0.42 


0.074 


FALSE 


FALSE 


Dwarf 


0.12 


0.002 


FALSE 


FALSE 


Dwarf 


2.21 


0.000 


FALSE 


FALSE 


Dwarf 


11.67 


0.000 


FALSE 


FALSE 


Dwarf 


8.41 


0.032 


FALSE 


FALSE 


Dwarf 


2.66 


0.150 


FALSE 


FALSE 


Dwarf 


0.60 


0.000 


FALSE 


FALSE 


Dwarf 


7.84 


1.000 


FALSE 


FALSE 


Dwarf 


14.02 


0.002 


FALSE 


FALSE 


Dwarf 


10.39 


0.042 


FALSE 


FALSE 


Dwarf 


3.35 


0.005 


FALSE 


FALSE 


Giant 


0.13 


0.010 


FALSE 


FALSE 


Giant 


90.38 


0.000 


FALSE 


FALSE 


Dwarf 


0.86 


0.000 


FALSE 


FALSE 


Dwarf 


4.67 


0.015 


FALSE 


FALSE 


Dwarf 


32.53 


0.906 


FALSE 


FALSE 


Dwarf 


59.89 


0.017 


FALSE 


FALSE 



the stars were identified and removed using tiie SysRem al- 
goritlim of lTamuz. Mazeii. fc Zuckeij (|2005l ). Tlie searcii for 
transit-like events and selection of astrophysically-plausible 
transit candidates was perfor med using the methodology of 
ICollier Cameron et al.l l|2006l ). backed up by visual inspec- 
tion of the ligh t curves to verify the reality of the transits 
as described bv [Christian et all (|2006h . 



In total these procedures yielded a "long-list" of 109 
transit candidates with periods less than 5 days, from which 
obvious impostors such as ellipsoidal variables and objects 
showing clear secondary eclipses on the best-fitting period 
had been eliminated. 
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Table 1 - continued 



ISWASP Tel/Inst Spec. Phot. RPM MS Prior Probability Select? Select? 

Class Class Class 5 = -2 In (Rp < l.bRj^p) (MCMC (MCMC 

V{M*,R,) only) & phot) 



Phot. Binaries: 








J003039.21+205719.1 


T,W 


SB2 


SE2 


J015711. 29+303447.7 


T 


SB2 


SE2 


J160242.43+290850.1 


W 




SE2 


J165423.72+241335.7 


W 




SE2 


J183104.01+323942.7 


S 


SB2 


SE2 


J183805.57+423432.3 


w 




SE2 


J205308.03+192152.7 


s 


SB2 


SE2 


J210909.05+184950.9 


w 




EV2 


J222353.83+412813.5 


w 




EV2 


J223651. 20+221000.8 


s 


SB2 


SE2 



Dwarf 


1.60 


0.022 


FALSE 


FALSE 


Dwarf 


1.42 


0.927 


TRUE 


FALSE 


Giant 


16.32 


0.000 


FALSE 


FALSE 


Giant 


2.88 


0.847 


TRUE 


FALSE 


Dwarf 


0.68 


0.997 


TRUE 


FALSE 


Giant 


4.08 


0.995 


TRUE 


FALSE 


Dwarf 


3.22 


1.000 


TRUE 


FALSE 


Dwarf 


0.36 


0.991 


TRUE 


FALSE 


Dwarf 


53.32 


0.000 


FALSE 


FALSE 


Giant 


1.13 


0.021 


FALSE 


FALSE 



Aperture blends: 








J025500.31+281134.0 


S,T 


< 10 m s-i 


AB 


J133156.81+460026.6 


W 




AB 


J161644.68+200806.8 


w 




AB 


J181454.99+391146.0 


w 




AB 


J184303.62+462656.4 


S 


I/S Nal, ClassIII 


AB 


J204142.49+075051.5 


w 




AB 


J204712.42+202544.5 


W,T 


< 1 km s~^ 


AB 


J204745.08+103347.9 


W 




AB 


J213416.37+205644.4 


W,T 




AB 


J215226.17+331424.7 


W 




AB 


J222317.60+130125.8 


s 


I/S Nal, ClassIII 


AB 


J223809.90+401038.1 


w 




AB 


J224104.57+363648.3 


w 




AB 



Giant 


15.19 


0.000 


FALSE 


FALSE 


Giant 


82.69 


0.000 


FALSE 


FALSE 


Giant 


79.66 


0.000 


FALSE 


FALSE 


Giant 


2.55 


1.000 


TRUE 


FALSE 


Giant 


60.26 


0.971 


FALSE 


FALSE 


Giant 


1.54 


1.000 


TRUE 


FALSE 


Giant 


1.93 


1.000 


TRUE 


FALSE 


Giant 


65.54 


0.780 


FALSE 


FALSE 


Dwarf 


35.47 


0.000 


FALSE 


FALSE 


Dwarf 


52.29 


0.036 


FALSE 


FALSE 


Giant 


38.09 


0.994 


FALSE 


FALSE 


Giant 


5.55 


1.000 


FALSE 


FALSE 


Giant 


43.77 


0.994 


FALSE 


FALSE 



Other blends, giants: 



J161732.90+242119.0 


S,E 


< 10 m s^i 




Dwarf 


0.53 


1.000 


TRUE 


TRUE 


J210318.01+080117.8 


S 


I/S Nal, BisVar 




Dwarf 


0.26 


0.983 


TRUE 


TRUE 


J005225.90+203451.2 


S,T 


< 10 m s~^ 




Dwarf 


14.03 


0.008 


FALSE 


FALSE 


J181252.03+461851.6 


S,E 


< 10 m s~^ 




Dwarf 


2.46 


0.036 


FALSE 


FALSE 


J204125.28+163911.8 


S 


ClassIII 




Giant 


37.85 


1.000 


FALSE 


FALSE 


J205027.33+064022.9 


S 


< 10 m 




Dwarf 


33.23 


0.986 


FALSE 


FALSE 


J214151.03+260158.5 


S 






Giant 


26.81 


1.000 


FALSE 


FALSE 



2.2 OHP 1.93-m spectroscopy 

The initial spectroscopic followup of transit candidates iden- 
tified from the 2004 SuperWASP survey was carried out 
during three 4- night observing runs on the 1.93-m telescope 
at I'Observatoire de Haute-Provence (OHP). The first two 
runs, in 2006 May and June, utilise d the ELODIE spec - 
trograph and radial-velocity software (jBouchv et al ] l2006t ). 
ELODIE was replaced by its successor SOPHIE in time for 
the third run of the series, in 2006 August. Both instruments 
were used in their "high-efficiency" modes, which cover the 
wavelength range from 387 to 694 nm with resolving power 
A/AA = 35000. Radial velocities are computed at the tele- 
scope, using cross-correlation with a mask spectrum, using 
ThAr arc spectra to monitor the stability of the instrument. 
The typical radial-velocity precisions were of order 40 m 
s~^ for narrow-lined objects observed with ELODIE, and 
12 m s~^ with SOPHIE. The cross-correlation functions al- 
low several types of astrophysical false positive to be iden- 



tified from a single exposure. Targets were rejected if the 
first-epoch spectrum revealed them to be double-lined or 
rapidly rotating (and hence presumably tidally locked by a 
stellar-mass companion). Examples of such cross-correl ation 
functions are illustrated in Fig. 1 of IPont et all (|2005al ). Ve- 
locity shifts of more than a few hundred m s~^ between the 
first spectrum and a second measurement on a subsequent 
night also led to rejection as a stellar binary. 

As Table [1] shows, 14 targets were eliminated as broad, 
single-lined or double-lined binaries during the ELODIE 
runs. 12 stars were identified as binaries with SOPHIE, while 
a further 10 were found to be either giants (deduced visu- 
ally from absence of pressure-broadening wings on the Na I 
D and Mg I b lines, and/or from the presence of narrow fore- 
ground interstellar sodium absorption), or to exhibit no de- 
tectable radial velocity variation (suggesting blending with 
a nearby eclipsing binary). 
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-lu -02 an u 



Figure 1. Blended eclipse of the planet candidate ISWASP 
J213416. 37+205644. 4. The top panel shows the low amplitude 
eclipse induced in the lightcurve by the eclipsing Algol BQ Peg 
at a distance of 69 arcsec. The bottom panel shows the lightcurve 
for BQ Peg itself. The lack of a sharp minimum on the primary 
eclipse is due to the poorer data quality as the flux approaches 
the SuperWASP faint limit at the V = 16.7 primary minimum. 
Both are phase folded on the BQ Peg period of 1.574520 days. 



tSWASP Jm i<iam2P56«*4 ffbnetQndtdatti 




UCOJ. I I . I I : , 

4« ; -0^ -u nil ibj I}* pid 

Figure 2. Aperture-blend test for planet candidate ISWASP 
J213416. 37+205644. 4. The plot shows the magnitude in the main 
photometric aperture minus that in the larger-outer aperture, 
phase folded on the determined period of the planet candidate. A 
V-shaped eclipse is clearly seen indicating this object is blended 
with a resolved eclipsing binary system, in this case BQ Peg. 



2.3 Tautenburg spectroscopy 

Further spectra of several targets were secured using 
the coude echelle spectrograph on the 2-m Alfred- Jensch- 
Teleskop at the Thiiringer Landessternwarte Tautenburg 
(TLS) over the winter of 2006-7. A 2-arcsec slit width was 
used, yielding a spectral resolving power A/ A A — 35000. 
For most observations, velocities were obtained by cross- 
correlation referenced to the ThAr wavelength calibration, 
using telluric lines as a secondary reference to monitor and 
correct possible instrumental shifts. The overall velocity pre- 
cision was typically of order 200 m s~^. The rejection criteria 
were similar to those described above for the OHP observa- 
tions. Six additional stars were eliminated as binaries on the 
basis of these observations. 



2.4 INT spectroscopy 

Spectra were secured of several previously-unobserved tar- 
gets using the Intermediate-Dispersion Spectrograph (IDS) 
on the 2.5-m Isaac Newton Telescope (INT) at the Obser- 
vatorio del Roque de los Muchachos, La Palma, during 2007 
April. Similar rejection criteria were applied. Five stars were 
found to be single-lined binaries. 

3 FURTHER PHOTOMETRIC TESTS 
3.1 Resolved aperture blends 

Eclipsing binary systems blended with nearby bright, non- 
variable stars represent the most common typ e of mimic in 
planetary transit surveys ijPont et al.l l2005ar i. The WASP 
data reduction pipeline performs photometry in three con- 
centric circular apertures centred on each object, with radii 



of 2.5, 3.5 and 4.5 13.5-arcsec pixels. In cases where an eclips- 
ing binary system lies just outside the main (3.5-pixel ra- 
dius) photometric aperture of a non-variable star, the two 
systems will contribute different amounts of light to differ- 
ent apertures. Light leaking from the wings of the binary's 
image into the apertures of the non-variable star may mimic 
a shallow transit at times when the binary is in eclipse. The 
amount of leakage is greatest in the outer apertures. As a 
consequence, the depth of the "transit" in the non-variable 
object is greatest when measured in the outer apertures. 

The eclipsing binary BQ Pegasi is a typical example. 
This is an Algol system with primary eclipses 2.6 mag deep, 
located 69 arcsec away from a 12th magnitude non-variable 
star (ISWASP J213416.37+205644.4). The centre-to-centre 
separation of these two systems is approximately 5 pixels 
causing a small amount of overlap between the photometric 
apertures. This overlap was sufficient to induce a 2% dip 
in brightness in the non-variable star, closely mimicking the 
photometric signal of a transiting planet with a period of 
1.57455 days - identical to BQ Peg. As a result the non- 
varying star was flagged by the transit detection algorithm 
as a high-priority planet candidate. The planet candidate 
lightcurve is shown in Figure [1] together with the lightcurve 
for BQ Peg. We see that the phase of the false transit, folded 
on the period of BQ Peg, matches that of the binary's pri- 
mary eclipse. Without knowledge of the variable nature of 
the blending system this candidate would be difficult to elim- 
inate without further foUowup. However, a plot of the differ- 
ence in magnitudes between the main and outer photometric 
apertures (Figure [2}, phase folded on the orbital period de- 
termined for the planet candidate, shows a clear V-shaped 
eclipse indicating that the star is blended. A non-blended 
system would be not be expected to produce a change in 
the flux ratio between different apertures in this situation. 

Thirteen objects failed the photometric aperture-blend 
test. Five among these were also observed spectroscopically. 
Two were found to be giants; of the two that were observed 
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more than once, neither showed significant radial-velocity 
variability. 



3.2 Ellipsoidal variables and secondary eclipses 

Because the transit search was restricted to periods less 
than 5 days, a few objects were detected in which secondary 
eclipses and/or ellipsoidal variability were present, but in 
which the best-fitting period found by the transit-search 
algorithm was half the true period. To guard against this 
possibility we examined the light curves of all candidates 
phase-folded on twice the best-fitting period. 

Two were eliminated photometrically as exhibiting el- 
lipsoidal variability on twice the orbital period found by 
the transit detection software. A further seven objects were 
found to exhibit unequal depths or displacements in phase 
between odd-numbered and even-numbered transits, indi- 
cating that the periods were indeed twice those found by 
the transit-search software, and that secondary eclipses were 
present. The five among these that were also observed spec- 
troscopically were all found to be double-lined binaries. 



4 MARKOV-CHAIN MONTE-CARLO 
MODELLING 

Markov-chain Monte-Carlo (MCMC) methods are rapidly 
gaining popularity as an efficient and informative means 
of solving of multivariate parameter-fitting problems in as- 
tronomy and many other branches of science. They pro- 
vide not just a means of optimising the fit of a model to 
data, but their mode of operation allows a precise explo- 
ration of the join t posterior probability distribution of the 
fitted parameters. iTegmark et al.l(|2004l ) give a clear descrip- 
tion of the technique and its use for deriving limits on cos- 
mological parameters from the spatial power spectrum of 
the co smic microwave background. iFord ()2006. 'l and .Gregorvl 
have applied MCMC to the derivation of orbital pa- 
rameters and d etection of ad d itiona l planets from radial- 
velocity curves. iHolman et al.l ()2006l ') have shown Markov- 
Chain Monte-Carlo (MCMC) methods also to be particu- 
larly well-suited to the problem of deriving the physical pa- 
rameters of star-planet systems by optimising model fits to 
the light curves of transiting exoplanets. 

In our implementation of MCMC, we assume the 
planet's orbit to be circular. We characterise the system 
using the six parameters {To, P, A_F, iT, &, M,}. Here To is 
the epoch of mid-transit, P is the orbital period, A_F is the 
fractional fiux deficit that would be observed during transit 
in the absence of limb darkening, It is the total duration 
of the transit from first to last contact, b is the impact pa- 
rameter of the planet's path across the stellar disc in units 
of the primary's radius, and M* is the mass of the primary 
in solar units. These six quantities constitute the "proposal 
parameters", which are allowed to perform a random walk 
through parameter space, generating a cloud of points that 
map out the joint posterior probability distribution (Fig. [H]). 
The best-fitting values of the proposal parameters are listed 
for all 67 stars in Table [2] 

At each step in the MCMC procedure, each proposal 
parameter is perturbed from its previous value by a small 



random amount: 



To,!. 


~ To.i- 


i+aroG(0,l)/ 




= P-i 


+ apG{0,l)f 


AF, 


= AF,- 


-1 +o-afG(0,1)/ 


tT,i 




i+at^G{0,l)f 


h 




i+atG{0,l)f 






+ (7A/G(0, 1)/ 



where G(0, 1) is a random Gaussian deviate with mean zero 
and unit standard deviation. The scale factor / is an adap- 
tive step-size controller whose value is initially set to 0.5, 
but which evolves together with the estimated parameter 
uncertainties as described in Section [4.11 below. 

The first four parameters (To, P, AF and tr) are more 
or less directly measurable from the folded light curve. Their 
initial values and their associated one-sigma uncertainties 
are taken directly from the results of the accelerated Bo x 
Least- Squares algorithm of ICollier Cameron et all l|200d ). 
The impact parameter is given an initial value bo = 0.5 
and a one-sigma uncertainty ab = 0.05. The stellar mass 
Mt is initially set to the value Mq derived from the J — H 
colour using the calibration described in Appendix [BJ and 
assigned an arbitrary but plausible one-sigma uncertainty 

GM = O.lMo. 

Once the physical parameters -R*, Rp, a and cosz have 
been derived from the proposal parameters AF, tr, b and 
P using the relationships presented in Appendix El it is 
straightforward to compute the projected separation of cen- 
tres (in units of the primary radius) at any time tj of obser- 
vation 



sin^ + {bR:,/aY 



(1) 



The orbital phase angle at time tj is (j>3 ~ '^''^{tj — To)/P. 

We compute the fiux defi cit at all obser v ed or bital 
phases using the algorithm of iMandel fc Agol (|2002h for 
small planets with the 4-coefficient limb-darkening model of 
[Claret (2000). The four limb-darkening coefficients are in- 
terpolated from Claret's tables for the R band (whose mean 
wavelength approximates the effective wavelength of the un- 
filtered WASP instrumental bandpass) using the stellar ef- 
fective temperature determined from the J — H colour via 
equation IBll We adopted surface gravity loggr = 4.5, a mi- 
croturbulent velocity Vtmh ~ 2 km s"'^, and a metallicity 
[Fe/H] — 0.1 on the grounds that planet-bearing stars tend 
to have heavy-element abundances slightly above solar. Af- 
ter converting these fiux deficits to model magnitudes /ij 
relative to the fiux received from the system outside transit, 
we compute the zero-point offset from the observed magni- 
tudes m,: 



. Ej("ii-W)wj 
Am = „ : 



(2) 



where the observational errors aj define the inverse- variance 
weights Wj = 1/cr^. We thus obtain the fitting statistic for 
the set of model parameters pertaining to the ith step of the 
Markov chain: 



ii{To,P,AF,tT,b,M,) = J2 



{rrij - - Amf 



(3) 
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Table 2. Observable and physical parameters for the 67 stars in the sample, computed by optimising the posterior probability Q 
incorporating the main-sequence prior constraint on the primary mass and radius. 



ISWASP Transit Orbital Transit Transit Impact Stellar Primary Secondary 

Epoch To period P duration depth parameter Mass radius radius 

(HJD) (days) tr/P AF b AL/Mq R*/Rq i?2/-Rjup 



Planets: 

J002040.07+315923.7 3219.5240 2.52065 0.055 

J160211. 83+281010.4 3178.2780 3.94164 0.032 

J203054.12+062546.4 3216.7138 2.15202 0.034 



0.010 0.028 1.24 1.38 1.36 

0.019 0.402 0.99 0.97 1.30 

0.015 0.696 0.84 0.83 1.02 



Sp. Binaries: 



J130322.00+350525.4 


3167.4755 


2.67382 


0.055 


0, 


.014 


0, 


.038 


1.05 


1.28 


1.48 


J133339.36+494321.6 


3175.1566 


3.64942 


0.038 


0, 


.009 


0, 


.021 


1.33 


1.19 


1.08 


J134815.27+464811.0 


3169.5603 


1.05603 


0.056 


0, 


.009 


0, 


.807 


1.04 


1.04 


0.98 


J174118.30+383656.3 


3197.8605 


4.24481 


0.037 


0, 


.010 


0, 


.280 


1.28 


1.32 


1.28 


J175401. 58+322112. 6 


3194.0863 


1.94918 


0.050 


0, 


.012 


0, 


.537 


1.13 


1.10 


1.15 


J181858.42+103550.1 


3170.8197 


2.46530 


0.056 


0, 


.010 


0, 


.113 


1.41 


1.39 


1.34 


J184119.02+403008.4 


3188.2544 


3.73765 


0.040 


0, 


.013 


0, 


.058 


1.06 


1.15 


1.29 


J203704.92+191525.1 


3215.3230 


1.68008 


0.045 


0, 


.009 


0, 


.762 


1.12 


1.11 


1.02 


J210009.75+193107.1 


3194.0029 


3.05562 


0.041 


0, 


.007 


0, 


.516 


1.61 


1.40 


1.13 


J211608.42+163220.3 


3219.5782 


3.46829 


0.031 


0, 


.015 


0, 


.842 


1.23 


1.35 


1.60 


J212855.03+075753.5 


3220.9833 


4.68994 


0.023 


0, 


.030 


0, 


.757 


0.95 


0.92 


1.56 


J223320.44+370139.1 


3240.4446 


1.87478 


0.050 


0, 


.010 


0, 


.565 


1.29 


1.17 


1.11 


J234318.41+295556.5 


3245.1886 


4.24098 


0.030 


0, 


.021 


0, 


.008 


0.78 


0.85 


1.20 


J010151. 11+314254.7 


3230.3097 


3.65098 


0.050 


0, 


.014 


0, 


.010 


0.76 


1.28 


1.45 


J023445.65+251244.0 


3234.8368 


1.55373 


0.082 


0, 


.018 


0, 


.020 


0.93 


1.26 


1.65 


J031103.19+211141.4 


3250.5698 


2.72925 


0.063 


0, 


.031 


0, 


.020 


0.94 


1.36 


2.35 


J051221. 34+300634.9 


3253.3499 


1.23758 


0.085 


0, 


.028 


0, 


.657 


1.03 


1.36 


2.22 


J115718.66+261906.1 


3156.3795 


2.45397 


0.058 


0, 


.014 


0, 


.066 


1.36 


1.38 


1.60 


J141558.71+400026.7 


3170.0586 


1.64277 


0.055 


0, 


.050 


0, 


.995 


1.33 


1.29 


3.95 


J152131.01+213521.3 


3154.5622 


4.01493 


0.043 


0, 


.026 


0, 


.082 


1.19 


1.32 


2.08 


J153135.51+305957.1 


3181.9327 


4.46754 


0.047 


0, 


.029 


0, 


.070 


0.94 


1.40 


2.31 


J165949.13+265346.1 


3200.5994 


2.68305 


0.042 


0, 


.023 


0, 


.885 


1.11 


1.48 


2.35 


J172549.13+502206.4 


3198.8106 


2.27129 


0.044 


0, 


.025 


0, 


.785 


1.01 


1.19 


1.84 


J173650.17+105557.9 


3193.1694 


3.44102 


0.041 


0, 


.030 


0, 


.037 


0.94 


1.03 


1.74 


J173748.98+471348.7 


3182.5662 


3.33791 


0.046 


0, 


.008 


0, 


.025 


0.86 


1.19 


1.05 


J174327.81+582512.7 


3213.7625 


2.84493 


0.064 


0, 


.013 


0, 


.016 


1.01 


1.53 


1.69 


J175511. 09+134731. 5 


3203.1171 


2.44367 


0.057 


0, 


.019 


0, 


.602 


1.03 


1.43 


1.92 


J175620.84+253625.7 


3204.0353 


4.41602 


0.034 


0, 


.031 


0, 


.667 


1.05 


1.28 


2.18 


J180010.55+214510.2 


3198.2836 


4.41513 


0.023 


0, 


.051 


0, 


.751 


0.67 


0.75 


1.66 


J182620.36+475902.8 


3205.5267 


3.04416 


0.086 


0, 


.068 


0, 


.935 


0.05 


0.91 


2.99 


J202824.02+192310.2 


3211.1395 


1.25788 


0.092 


0, 


.018 


0, 


.016 


1.30 


1.36 


1.76 


J203906.39+171345.9 


3209.4011 


2.69705 


0.045 


0, 


.021 


0, 


.758 


1.16 


1.42 


1.99 


J215802.14+253006.1 


3215.6417 


1.51272 


0.046 


0, 


.018 


0, 


.332 


0.75 


0.67 


0.87 


J231533.56+232637.5 


3225.8489 


4.56289 


0.061 


0, 


.016 


0, 


.003 


0.35 


1.39 


1.74 



4.1 Main-sequence prior 

In Bayesian terms, the likelihood of obtaining the observed 
data D given the model defined by a particular set of pro- 
posal parameters is 

P(D\To,P,AF,tT,b,M,)cxeM-xV^), (4) 

but the full posterior probability distribution for the data 
and the model depends on the prior probability distribution 
for each of the model parameters. We are only interested in 
solutions for which the companion is a planet-sized object 
yielding a transit of observable depth, so we restrict the im- 
pact parameter to the range < fo < 1.0, rejecting proposal 
steps that fall outside this range. 



For most of the remaining parameters the uniform prior 
implied by the random-walk nature of MCMC is valid, but 
the stellar mass and radius are already determined a pri- 
ori from the J — H colour, under the implicit assump- 
tion that the star is single and on the main sequence. Un- 
der this prior assumption we expect the stellar mass to 
lie somewhere within an approximately gaussian distribu- 
tion with mean AIq and standard deviation aM = A/o/10 
(i.e. the same arbitrary but plausible value used to deter- 
mine the average jump size in M,). We use a power-law 
approximation to the main-sequence mass-radius relation 
to define a prior probability distribution for i?* with mean 
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Table 2 - continued 



IsWAsP 


Transit 


Orbital 


Transit 


Transit 


Impact 


Stellar 


Primary 


Secondary 




Epoch To 


period P 


duration 


depth 


parameter 


Mass 


radius 


radius 




(HJD) 


(days) 


tr/P 


AF 


b 


M^/Mq 


R*/Rq 




Phot. Binaries: 


















J003039. 21+205719.1 


3230.6250 


4.56617 


0.028 


0.020 


0.792 


1.21 


1.34 


1.86 


J015711. 29+303447.7 


3237.7106 


2.04321 


0.056 


0.012 


0.428 


1.05 


1.17 


1.27 


J160242. 43+290850.1 


3187.2622 


1.30433 


0.082 


0.044 


0.681 


0.66 


1.13 


2.30 


J165423. 72+241335.7 


3188.8269 


2.57122 


0.045 


0.031 


0.030 


0.60 


0.80 


1.36 


J183104.01+323942.7 


3210.9015 


2.37843 


0.040 


0.008 


0.769 


1.30 


1.33 


1.18 


J lS3o05. 57+423432. 3 


3199.0823 


3.51709 


0.036 


0.016 


0.046 


0.59 


0.82 


0.99 


J 205308. 03+192152. 7 


3204.4935 


1.67630 


0.076 


0.005 


0.008 


1.22 


1.41 


1.00 


J 210909. 05+184950. 9 


3216.9989 


2.91766 


0.044 


0.007 


0.470 


1.70 


1.45 


1.16 


J 222353. 83+412813. 5 


3239.3999 


1.55304 


0.130 


0.038 


0.112 


0.48 


1.51 


2.87 


J 223651. 20+221000. 8 


3205.6264 


3.22659 


0.032 


0.056 


0.865 


0.77 


0.91 


2.42 


Aperture blends: 


















J025500.31+281134.0 


3237.4163 


2.16594 


0.068 


0.029 


0.034 


0.70 


1.14 


1.88 


J133156. 81+460026.6 


3166.2329 


3.16678 


0.091 


0.063 


0.072 


0.08 


0.91 


2.22 


J161644.68+200806.8 


3243.4648 


3.96830 


0.069 


0.178 


0.994 


0.10 


0.85 


4.90 


J181454.99+391146.0 


3192.6921 


1.10269 


0.067 


0.018 


0.081 


0.47 


0.65 


0.85 


J184303. 62+462656.4 


3192.1856 


3.33871 


0.060 


0.020 


0.066 


0.22 


0.94 


1.30 


J204142. 49+075051. 5 


3207.5200 


2.76243 


0.037 


0.008 


0.110 


0.59 


0.75 


0.65 


J204712. 42+202544.5 


3208.9201 


2.61303 


0.041 


0.020 


0.004 


0.60 


0.76 


1.06 


J204745. 08+103347.9 


3200.9224 


3.23610 


0.062 


0.021 


0.061 


0.18 


0.88 


1.23 


J2iv)4lD.o7+zu5u44.4 


3220.5575 


1.57455 


0.105 


0.019 


0.013 


0.63 


1.42 


1.89 


J 215226. 17+331424. 7 


3226.5888 


1.06591 


0.149 


0.018 


0.039 


0.41 


1.34 


1.75 


JZZZolY.DU+ioUlzO.O 


ozz7.4z5U 


1.61385 


0.085 


0.016 


0.060 


0.39 


1.00 


1.23 


J 223809. 90+401038.1 


3240.3573 


1.40652 


0.065 


0.016 


0.023 


0.52 


0.77 


0.94 


J 224104. 57+363648. 3 


3245.1212 


1 '70'irio 

1.73399 


0.081 


0.019 


0.100 


0.31 


0.92 


1.24 


Other blends, giants: 


















J161732. 90+242119.0 


3182.9227 


1.45378 


0.055 


0.012 


0.436 


0.72 


0.82 


0.87 


J210318. 01+080117.8 


3214.3962 


1.22395 


0.077 


0.013 


0.255 


1.06 


1.10 


1.22 


J005225. 90+203451. 2 


3230.5589 


1.71894 


0.059 


0.031 


0.869 


0.81 


1.27 


2.36 


J181252. 03+461851. 6 


3210.2889 


2.52501 


0.039 


0.033 


0.831 


0.91 


1.10 


1.98 


J204125. 28+163911. 8 


3214.3614 


1.22154 


0.096 


0.006 


0.027 


0.38 


0.97 


0.74 


J205027.33+064022.9 


3209.8236 


1.22927 


0.125 


0.007 


0.007 


0.74 


1.58 


1.31 


J214151. 03+260158.5 


3238.9321 


1.82588 


0.064 


0.007 


0.031 


0.38 


0.84 


0.69 



Hastings rule; if Qi < Qi-i the new set is accepted; if on 
the other hand Qi > Qi-i, the new set is accepted with 
probability exp(— A(3/2), where AQ = Qi — Qi-i- The algo- 
rithm first converges to, then explores the parameter space 
around a constrained optimum solution that represents a 
compromise between fitting the hght curve and reconciling 
the resulting stellar dimensions with prior expectations de- 
rived from the J — H colour. As we shall see later, the value 
of the main-sequence prior ViM^^i, R,,i) at the global mini- 
mum of Q constitutes a useful measure of how far the stellar 
parameters must be displaced from the main sequence val- 
ues appropriate to the star's colour in order to fit the transit 
light curve. 

We introduce a step-size controller / which is adjusted 
on every 100th step to ensure that t he acceptance rate is 
held close to the optimal value of 0.25 ijTegmark et al.ir2004] l. 
The value of / is determined by the simple linear algorithm 
fnew = 400/oid/j, where j is the total number of propos- 



.Ro = Mf? * (|Tinglev fc Sacketd lioosi ) and hence standard 
deviation an = 0.8{Ro /Mo)crM ■ 

This gives a joint prior probability distribution for the 
values of the proposal parameterM* and the derived physical 
parameter i?« of the form 

W„i?„0=exp ^^)-^'^ 

Since the posterior probability distribution is 
V{M,,i,R,,i)V{D\To,P,AF,tT,b,M,), we impose the 
prior on Mt_i and by replacing with the logarithm 
of the posterior probability distribution 

Q,.X-+ ^^*%^°^% ^^^^ (6) 

as the statistic on which acceptance of a set of proposal 
parameters is decided. 

For every new proposal set generated, the decision as to 
whether or not to accept the set is made via the Metropolis- 
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Figure 3. The lower triangular matrix shows the correlation diagram for the six proposal parameters {Tq, P, AF, ty, 6, M*} for an 
MCMC analysis of the 2004 SuperWASP light curve of the planet host star XO-1. The MCMC run imposed Bayesian priors on both the 
stellar mass and the stellar radius, derived from the J — H colour and the main-sequence mass-radius relation. The shape of the transit 
profile allows a wide range of intermediate values for the impact parameter b. Nonetheless, the majority of the proposal parameters are 
very close to being mutually uncorrelated, with the exception of b and AF, which are correlated through the effect of limb darkening. 
The upper triangular matrix shows the mutual correlations among the physical parameters b and Rt/a and Rp / ij« . The well-known 
correlation between b and Rt/a is very much stronger than those seen among the six observational parameters that we have adopted 
here. 



als generated and tested during the previous 100 successful 
steps. We find that values in the range 0.5 < / < 1 achieve 
the necessary acceptance rate when the correct values of 
the parameter uncertainties are used. In any case, we find 
that all six parameters invariably converge to their optimal 
values within 500 steps. Once this "burn-in" phase is com- 
plete, the uncertainties on the data are adjusted to take 
account of correlated errors by rescaling to ensure that the 
value of the optimal solution is equal to the number 

of degrees of freedom. The parameter uncertainties are sub- 
sequently allowed to evolve as the computation progresses, 
being recomputed from the Markov chains themselves ev- 
ery 100 successful steps. This ensures that the step size in 
each dimension samples the parameter space adequately. Af- 
ter discarding the first 500 steps we compared the variances 
within and bet ween five independent su b-chains, using the 
test statistic of iGelman fc RubinI (|1992| ). The value of the 
Gelman-Rubin statistic was invariably within a fraction of 
one percent of unity, verifying that the chains were properly 
converged and well-mixed. 

For all its advantages, MCMC is a discrete, stochastic 
process. After each MCMC run was completed, we there- 
fore refined the values of the six proposal parameters that 
represent the optimum solutio n for Q using the AMOEBA 
(downhill simplex) algorithm l|Press et al.lll993 ), using the 



best-fitting step of of the Markov chain as one of the ini- 
tial vertices and the first six successful steps of the chain 
as the others. The derived primary and secondary radii (in 
solar and jovian units), are listed alongside the best-fitting 
proposal parameters for all 67 stars in Table 



4.2 Orthogonality properties 

We note that the proposal parameter set 
{To,P,AF,tT,b,M,} is particularly well-suited to MCMC, 
because the first four parameters have approximately 
Gaussian probability density functions that are directly 
related to observationally-determined quantities. As Fig. [3] 
illustrates, all six are close to being mutually uncorrelated. 
This circumvents the difficulty commonly encountered 
in some applications of MCMC to transit modelling, 
in which the normalised stellar radius R^/a is used di- 
rectly in place of tr- The parameters Rp/Rt, R^/a and 
b are strongly covariant, leading to correlation lengths 
of several hundred t o 1000 steps in the Markov chain 
iHolman et all l2006l : IWinn. Holman. fc Roussanoval |2007| ) 
if they are used directly as proposal parameters. Shorter 
correlation lengths are desirable, since they allow more 
statistically-independent samples of the parameter space 
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Figure 4. Impact parameter versus planet radius and stellar radius versus stellar mass for MCMC analysis of the SuperWASP light 
curves of the grazing broad-lined spectroscopic binary ISWASP J215802. 14+253006.1. The left-hand realisation imposes a prior on the 
stellar mass, derived from the J — H colour. The right-hand panels show the effect of imposing an additional main-sequence prior on the 
stellar radius. 



to be generated in less computing time. In our parameter 
space the equivalent observational parameters AF, It and 
h are only very weakly covariant. Their use reduces the 
correlation lengths in the chains for these parameters to 
between 4 and 20 steps, yet they yield Rf/a and R^/Ri^ 
direct ly via eqs. lAll and IA2I We note that iBurke et al.l 
have arrived at a very similar conclusion, and that 
further shortening of the correlation length is possible if 
their orthogonalisation procedure is employed. 



5 SECONDARY RADII AND IMPACT 
PARAMETERS 

The use of the 2MASS J — H colour to place prior con- 
straints on the stellar mass and radius is of key importance 
in extracting essential information about the system dimen- 
sions from SuperWASP hght curves. The value of the prior 
at the constrained optimum solution gives a powerful diag- 
nostic for the presence of multiple stellar spectra, and for 
the host star's proximity to the main sequence. 

To illustrate this point we show in Figs and [S] 
the MCMC realisations of the relationship between im- 
pact parameter and planet radius for WASP-1, WASP-2 
and the broad double-lined spectroscopic binary ISWASP 
J215802. 14+253006.1. We performed two MCMC analyses. 
The first imposed a prior on the stellar mass only. The poste- 
rior probability density is then proportional to exp(— Qi/2), 



where 

_ 2 , {M„i-Mof 

- X. + 

is used as the statistic to which the Metropolis-Hastings 
rule is applied. The second analysis imposed the full main- 
sequence prior as described in Eq. [6] 

ISWASP J215802. 14+253006.1 was selected as a candi- 
date, having been found to exhibit transits 0.024 magnitudes 
deep. Its 2MASS J — H = 0.29 suggests an early G spec- 
tral type. It passed the t ransit depth and duration tests of 
iTinglev fc SackettI (|2005h and all other photometric selec- 
tion tests. A single SOPHIE spectrum, however, revealed 
the presence of two rotationally-broadened stellar spectra. 

The left-hand panel of Fig. H shows the MCMC joint 
posterior probability distributions of the impact parame- 
ter and secondary radius, and of the primary mass and ra- 
dius of ISWASP J215802. 14+253006.1. The MCMC results 
for which the prior is applied to the the stellar mass only 
show that the impact parameter is high, indicating a graz- 
ing eclipse. The mass-radius plot reveals a further inconsis- 
tency: the inferred primary radius is twice that expected for 
a main-sequence star. Imposing the further Bayesian con- 
straint on the primary radius does little to improve matters. 
The prior's attempt to pull the primary's dimensions toward 
the main sequence is thwarted by the need to fit the V shape 
of the transit light curve. The impact parameter is reduced 
somewhat, but the primary mass is pulled toward a value 
too low to be consistent with the colour. The primary's ra- 
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Figure 5. As for Fig.|4] for the known exoplanet host star WASP-1. 
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Figure 6. As for Fig. 2] for the known exoplanet host star WASP-2. 
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Table 3. Comparison of stellar and planetary dimensions derived 
by other authors from high-precision photometry with those de- 
rived from SuperWASP photometry in the present study using 
the constrained MCMC analysis. The middle columns show the 

radius values of Shporer et al and Charbonneau et al. rescaled as 
1/3 

MJ to account for the differences in the other authors' assumed 
and our fitted stellar masses. 



Parameter Shporer et al Charbonneau et al This study 
(2007) rescaled (2007) rescaled 



J — H is c onsistent with the range of values determined by 
IStempel"s~et al. (2007 ) from their recent spectroscopic abun- 
dance analysis of WASP-1. The inflated radius found by 
Shporer et al. and Charbonneau et al. is consistent with 
the rather low log g — 4.28 ± 0.15 found by Stempels et al. 
WASP-1 thus provides valuable verification that the selec- 
tion criteria suggested in Section |S] below are not so strict as 
inadvertently to exclude inflated planets orbiting slightly- 
evolved stars. 



WASP-1: 
R* / Rq 

Rp/Rjup 
b 

M./M0 
WASP-2: 

R,/Rq 

Rp/Rjup 

b 

M./M0 



1.46 ± 0.06 
1.44 ±0.06 
0.03 ± 0.17 
1.24 (matched) 



1.490 ±0.033 
1.480 ±0.040 
< 0.336 
1.24 (matched) 



0.830 ± 0.033 
1.059 ±0.051 
0.731 ±0.026 
0.84 (matched) 



1 QC9+0.047 
i-Jo^-O.iie 

1 ■?';8+°°''* 

i.OOO_g j^Q4 

215+0 °9l 
'^•^^^''-0.156 

1 240- 12 
^•^^-0.17 



U.001_Q ggg 

1 017+°"^ 
0.753lg:;^? 
0.840_:S^,, 



dius remains too large for its mass. The probability that the 
secondary has a radius less than 1.5 times that of Jupiter 
remains negligible, while the statistic representing the main- 
sequence prior is driven to the very high value 5* — 32.53. 
The star ISWASP J215802. 14+253006.1 thus fails two im- 
portant statistical tests, and can be eliminated as a viable 
exoplanet transit candidate without the need for spectro- 
scopic follow-up. 

For those systems that appear to have single main- 
sequence primaries, the mass-radius constraint also serves 
to reduce the uncertainty in the impact parameter of the 
transit. By imposing a prior that pulls the stellar radius to- 
ward the main-sequence mass-radius relation, we reduce the 
degeneracy between the stellar radius and impact parameter 
that arises from the transit duration. Since the ratio of the 
secondary's radius to that of the primary depends on the 
transit depth and impact parameter, the uncertainty in the 
secondary's radius is substantially improved. 

In Table 13] we compare the stellar and planetary dimen- 
sions inferred from the SuperWASP li ght curves using this 
method, with the values determined by Shporer et al.l (|2007l ) 
from h igh-quality 7-band t ransit photometry of WASP-1, 
and bv lCharbonneau et al.l ()2007l ) from high-quality z-band 
observations of both WASP-1 and WASP-2. The imposition 
of the main-sequence prior fixes the impact parameter at 
values for both stars that agree well with those of both sets 
of authors. Our values for the radii of both planet and host 
star in the WASP-2 system are also in excellent agreement 
with those of Charbonneau et al. 

Our radii for WASP-1 and its planet are systemati- 
cally smaller than those found by Charbonneau et al and 
by Shporer et al. These authors both determine a pri- 
mary radius that is consistent with the star being some- 
what above the main sequence and hence slightly evolved, 
whereas our method artificially pulls the radius toward the 
expected main-sequence value. The mass that we derive from 



6 SELECTION OF EXOPLANETARY 
TRANSIT CANDIDATES 

Instruments such as SuperWASP are capable of detecting 
the transits of gas-giant planets orbiting stars on or at most 
slightly above the main sequence. The J — H colour of a 
viable candidate should thus yield a stellar mass and radius 
that lie near the main sequence, and which are consistent 
with the duration of the transit. The value of the main- 
sequence prior P(M, , 7?«) mecisures the displacement of the 
fitted stellar mass and radius at the global minimum of Q 
from the main sequence values we expect from the star's 
J — H colour. 

The radius derived for the companion should be con- 
sistent with its being a planet. Although the SuperWASP 
light curves are not always of sufficient quality to discern 
the shape of the transit, the main-sequence prior helps to 
break the well-known degeneracy between the impact pa- 
rameter b and the stellar radius i?* (which is strongly cor- 
related with the companion radius via the transit depth). 
The use of MCMC allows us to determine directly the prob- 
ability that the companion radius is less than some specified 
amount. 

For convenience we define the quantity 



S = -2\nT{M,,R,) 

as a measure of the discrepancy between the constrained 
optimum stellar dimensions and the values derived from J — 
H. 

In Fig. [7]we plot S against the probability that the com- 
panion has a radius less than 1.5 times that of Jupiter, for 
the three planet-host stars and the various classes of astro- 
physical false-positive in our survey sample. The diagram 
shows clearly that the sample is cleanly partitioned into ob- 
jects with high and low probabilities of having Jupiter-sized 
secondary components. 

There is a tight cluster of objects in the lower-right 
corner of Fig. [T] with the main-sequence prior giving S < 1 
and probabilities greater than 0.95 that the secondary radius 
is less than 1.5 times that of Jupiter. This cluster includes 
all three known planets in the sample, together with two 
unresolved aperture blends and ten spectroscopic binaries, 
two of which can be eliminated photometrically as exhibiting 
ellipsoidal variations or unequal odd- and even-numbered 
transits. 

In selecting targets for spectroscopic foUowup it would 
be unwise to set thresholds as tight as the bounds of this 
cluster. The light curve of the recently-discovered TrES-3 
(jO'Donovan et all [ioOTh is a good illustration of the need 
for caution. The light curve of TrES-3 was recorded by 
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Figure 7. Logarithm of main-sequence prior probability versus probability that companion radius is less than 1.5-Rjup for known planets 
(red circles) and other forms of astrophysical false positive. Open circles and diamonds represent spectroscopic binaries and unresolved 
aperture blends, which could only be identified spectroscopically as false positives. Upright and diagonal crosses are eclipsing binaries 
and resolved aperture blends, detected using the methodology of Section |3 . 21 and |3. II respectivelv. The right-hand panel gives an enlarged 
view of the lower-right corner of the left panel. 



the SuperWASP cameras as ISWASP J175207.01+373246.3 
and its transits were detected with the correct period 
ijLister et al. 2007) even though it did not satisfy our original 
selection criteria for foUowup at high priority. Applying the 
MCMC test to the SuperWASP light curve of TrES-3, we 
find that S — 0.06, indicating a stellar mass and radius con- 
sistent with expectations, but V{R2 < 1.5f?jup) — 0.236. 
The reason for the apparently low probability of the sec- 
ondary having a planet-like radius arises because TrES-3 
has a high impact parameter. Limb darkening produces a 
strong positive correlation between 6 and Rp, giving a range 
of solutions with 0.8 < & < 1 and 1.3 < Rp < 2.0i?jup. By 
requiring candidates for followup spectroscopy to have S < 5 
and V{R2 < 1.5i?jup) > 0.2, we leave a sufficient margin to 
allow for the possibility of detecting infiated planets tran- 
siting slightly evolved main-sequence stars at high impact 
parameters. 

Six of the 19 binaries that satisfy these more relaxed se- 
lection criteria are represented by diagonal crosses in Fig. [T] 
having been identified photometrically as mimics from their 
ellipsoidal variations or secondary eclipses. Similarly, three 
of the five surviving blends, denoted by upright crosses in 
Fig.O, fail the photometric aperture-blend test. They could 
thus have been eliminated without requiring spectroscopic 
follow-up. With the wisdom of hindsight, we would thus re- 
tain 18 candidates from our original sample of 67 stars for 
spectroscopic followup, of which two are unresolved aperture 
blends, thirteen are stellar binaries and three are planet host 
stars. The final selection decisions are given in the last col- 
umn of Table [1] 



7 DWARF-GIANT SEPARATION 

If proper motion information is available, a useful secondary 
check can be made on the luminosity class of the primary. 
Evolved giant and sub-giant stars have very similar photo- 
metric colours to their main-sequence counterparts and are 
frequently a source of contamination in exoplanetary transit 
surveys. They can easily be identified, however, from their 




Figure 8. Reduced proper motion against J — H for dwarfs (log 
g 4.0) and giants (log $C 3.0). 



proper motions. Giant stars are significantly more distant 
than dwarfs of the same magnitude, and so exhibit substan- 
t iall y smaller proper motion s . 

IGouM fc MorganI (|2003l ) found that a plot of Reduced 
Proper Motion (RPM) against Tycho-2 Bt - Vt pro- 
vided an effective means of separating giants from dwarfs. 
We adapted this method using the J — H colour index 
and p roper motions from the all-sky USNO-Bl.O cata- 
logue l|Monet et al.ll2003l ). Our model is calibrated using 
stars from the catalogu e s of Valenti fc Fischeij (|2005l ) and 
ICavrel de Strobel et all (|200l[ ) which list over 2000 FGK 
stars with high resolution, high signal-to-noise spectra from 
which the surface gravity has been determined to better than 
± O.ldex. The stars were cross-matched with the USNO- 
Bl.O catalogue to determine proper motions and with the 
2MASS catalogue to determine J—H using a simple nearest- 
neighbour search within a radius of 10 arcsec. The stars were 
partitioned into dwarfs with logg ^ 4.0 and giants with 
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Figure 9. SOPHIE spectra of the Mg I b triplet region in the 
K giant ISWASP J204125. 28+163911. 8 (upper) and the K dwarf 
planet-host star WASP-2. Note the difference in pressure broaden- 
ing of the wings of the Mg lb triplet lines at 5167. SA, 5172.7Aand 
5183.6A. 



logg ^ 3.0. These conservative limits reduce the chances 
of accidental elimination of mildly evolved dwarfs. This re- 
sulted in 1526 dwarfs and 1145 giants. Although this is a 
much smaller sample than that used by iGould fc Morgan! 
l|2003ll , the luminosity class determination is based on more 
accurate spectroscopic measurements rather than broad- 
band photometric calibrations. The Reduced Proper Mo- 
tion, Hj, is defined as: 

Hj = J + 51og(M) (7) 

where J is the 2MASS J apparent magnitude and /x is the 
proper motion. In Figure [5] we plot Hj against J — H. A 
clean separation is seen between the giants and dwarfs. As 
predicted the giants, in general, have lower Reduced Proper 
Motions than the dwarfs with the same J — H colour. A 4th 
order polynomial has been fitted to a subjectively-drawn 
boundary showing the best separation between the two lu- 
minosity classes where: 



y = -141.25(J-//)''-)-473.18(J-//)^ 

-583.6(J - Hf + 313.42(J - - 58 

The greatest cross-contamination occurs for early K- 
dwarfs/giants {J — H ^ 0.5) in agreement with Gould and 
Morgan's findings. Unfortunately, this sample lacks large 
numbers of K and M dwarfs due to the inherent difficulty in 
obtaining high signal-to-noise spectra for these faint stars. 
As a result, the dwarf/giant determination becomes increas- 
ingly uncertain with decreasing temperature. 

We stress that although the reduced proper motion test 
is a very useful secondary check on the MCMC statistics, 
we do not use it as a primary selection criterion. It carries 
too high a risk of inadvertently rejecting genuine K dwarfs 
which just happen to have low proper motions. Nonetheless, 
we note that all the objects classified as giants on the ba- 
sis of their reduced proper motions fail one or other of the 
selection criteria. In particular, we note from Table [1] that 
the photometric aperture blends are almost invariably also 
classified as giants. The RPM giants among the aperture 
that were observed with SOPHIE showed the low pressure- 
broadening in the line wings expected of giants (Fig. [9} . 



8 DISCUSSION AND CONCLUSIONS 

We conclude that MCMC analyses of transit light curves 
can be carried out in such a way as to break the degeneracy 
between impact parameter and stellar radius, if a Bayesian 
prior is imposed to constrain the stellar mass and radius to 
values consistent with the host star's J — H colour. The pos- 
terior probability distributions for the stellar and planetary 
parameters are obtained directly from the Markov chain. 
They provide a straightforward estimate of the probability 
that the secondary has a radius consistent with planetary 
status. 

The prior carries the implicit assumption that the star 
is single and on the main sequence. We have found in this 
study that in a large fraction of astrophysical false positives 
for which these assumptions are invalid, the value of the 
prior at the optimum solution (in the sense of minimising 
the quantity Qi defined in Eq. [HJ is significantly greater 
than unity, if the uncertainty in the stellar mass derived from 
J—H is fixed (arbitrarily but reasonably) at 10 percent. The 
three known planets in our sample and the newly-discovered 
TrES-3 all yield values of the prior less than unity at their 
optimal solutions. 

The probability that the secondary is of planetary di- 
mensions and the value of the prior at the optimum so- 
lution provide the two-dimensional classification scheme 
shown in Fig. [T] For the sample of stars studied here, which 
had already passed the prelim inary selection procedures of 
ICollier Cameron et al.l ()2006l l. this classification scheme is 
effective at identifying 75 percent of all types of aperture 
blend. The figure improves to 90 percent if the flux ratio in 
different photometric apertures is used to identify blends of 
non-variable stars with eclipsing binaries near the edge of 
the photometric aperture. 

The Bayesian selection scheme also eliminates roughly 
67 percent of all spectroscopic binaries. It is particularly ef- 
fective at removing grazing stellar binaries, most of which 
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yield very low probabilities of having planet-sized secon- 
daries, and many of which also produce values of the main- 
sequence prior that are too high to be consistent with a 
single, main-sequence primary of the observed J — H colour. 
Grazing binaries frequently show either ellipsoidal variations 
or unequal eclipses when folded on twice the best period 
found by the accelerated box least-squares algorithm, giv- 
ing an effective second line of defence against this type of 
impostor. 

The remaining candidates have demonstrably low com- 
panion radii, generally low impact parameters, and pri- 
maries that appear to be close to the part of the main se- 
quence suggested by their observed J — H colours. Many 
of them are likely to be astrophysically interesting binaries, 
whose Jupiter-sized stellar or substellar companions will add 
to our knowledge of the mass-radius relation at and below 
the bottom of the main sequence once their spectroscopic 
orbits are established. Other teams pursuing similar stud- 
ies h ave already achie ved impor tant advances in this area 
ifBouc hv ct al. 2005; Pont et all i2005a .b, 2006,). In planet- 
hunting terms alone, however, the prospects for a high "hit 
rate" in future follow-up campaigns are good: in our sample, 
confirmed transiting planets comprise nearly 20 percent of 
the survivors. 
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Figure 10. Geometry of partial transit phases, showing the angles a and /3 used in computing the fraction of the stellar disc obscured 
by the companion. 
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APPENDIX A: OBSERVABLE AND PHYSICAL 
PARAMETERS 



ISeager fc Mallen-Ornela^ ()2003l l characterised the light 
curve of a transit candidate in terms of four observable quan- 
tities: the transit depth A_F expressed as a fraction of the 
system flux outside transit, the total transit duration tr 
from first to fourth contact, the duration tp of the flat part 
of the transit from second to third contact, and the orbital 
period P. Assuming that the transit is total, that all the light 
emanated from the primary star, and that the orbit was cir- 
cular, they derived expressions for the ratio of the planet's 
radius Rp to that of the star _R«, the ratio of R, to the or- 
bital separation a, and the impact parameter b = acosi/i?*. 
Kepler's third law then gives the stellar density p/pQ, under 
the further assumption that the mass ratio Mp/M^ << 1. 

Here we adopt a similar approach, generalised to 
allow for grazing as well as total eclipses. Following 
ISeager &: Mallen-Ornela^ l|2003h and using the approxima- 
tion that Rp -\- Rt « a we determine the stellar radius 
from the transit depth, duration and impact parameter via 
the relation: 



tr 



62 



(Al) 



For transits where the companion is fully silhouetted against 
the primary during the middle part of the transit we neglect 
limb darkening and deflne 



(A2) 



For partial transits, or at phases where the projected sepa- 
ration z of the centres lies in the range 1 ~ Rp/R, < z < 



1 + Rp/Rt, the fraction of the primary's disc that is obscured 
by the companion is 

R 

AF = — ^ (tt - /3 + sin /3 cos /3 + Q - sin Q cos a) . (A3) 
R* 

The angle a is subtended at the centre of the primary by 
the line of centres and the radius of the primary at the point 
where the two limbs intersect. The angle /3 is subtended at 
the centre of the planet by the line of centres and the radius 
of the planet at the point where the two limbs intersect 
(Fig, nop . The angles a and (3 are deflned by 



3 1 
cos a — -z -\ 

and 



{Rp/ Rt) 



2z 



cos /3 = 



{Rp/R* 



2zRp/R, 



(A4) 



(A5) 



For cases where Rp/ R, < 0.3 or so, the fractional obscu- 
ration of the primary's disc can be approximated to better 
than 4 percent by a linear function of the separation z of 
the centres: 



2R, 



1 + 

2R, 



(A6) 



This equation is quadratic in Rp/R, and has the positive 
solution 



(A7) 



Rp _ z-l + ^{l-zP+8AF 
R,~ 2 

Since z — b a,t mid-transit, the ratio of the planetary and 
stellar radii can be determined from the transit depth AF 
at mid-transit even for grazing eclipses, provided b is known. 

In practice, the combination of observational errors and 
the effects of limb darkening make it very difficult to deter- 
mine the duration tp oi the middle phase of the transit with 
any degree of reliability. For this reason we parametrise the 
form of the transit profile i n term s of AF, tr, b and P. Like 
ISeager fc Mallen-Ornela^ (|2003h we use Kepler's third law 
to determine the stellar density in solar units: 



= 0.0134063 ( ' 



R*J Vlday 



(A8) 
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In order to close the system, we need either an addi- 
tional equation or an independent determination of the stel- 
lar mass. It is then straightforward to determine the stellar 
radi us (in solar units) from th e stella r density. 

ISeager fc Mallen-Ornela^ (|2003l ) chose the former ap- 
proach, imposing a main-sequence mass-radius relation on 
the primary of the form i?« cx M^'^. We choose instead to 
use the J — H colo ur index taken fro m the 2MASS Point 
Source Catalogue o f ICutri et all (j2003l ) to estimate the stel- 
lar effective temperature and hence the corresponding main- 
sequence stellar mass, as described in Appendix iBl 



APPENDIX B: RADIUS ESTIMATION USING 

J-H 

The 2MASS J — H colour index was calibrated using stellar 
temperatu res from a sample of 1 00,000 Tycho-2 FGK dwarf 
stars from lAmmons et al.l (|2006t ). selected to assist in target 
selection for radial-velocity planet searches. The tempera- 
tures were calculated by fitting spline functions to broad- 
band Hipparcos/Tycho-2 and 2MASS photometry using a 
training set of stars observed with the Keck HIRES (High 
Resolution Echelle Spectrometer). The temperature model is 
quoted with an accuracy of better than lOOK. We extracted 
a subset of approximately 65000 stars from this sample for 
which the errors in the effective temperature were less than 
150K and the photometry errors less than 1 percent. The 
J — H index was then plotted against temperature (Fig. IBip 
and fitted with a linear relation 



Tcff = -4369.5( J - H)+ 7188.2, 



(Bl) 



which is valid over the approximate temperature range 
4000 < Teff < 7000 K. The J-H index is preferred over 
J — K as the 2MASS errors for H magnitude are generally 
lower than for K. 

The scatter of the data about the linear fit is 114K and 
is an improvement of the fit given by the relations of 
l|2000l) . 

The stellar radius is then calculated from a polynomial 
fit to the temperature/radius relat ion for mai n-sequence 
stars tabulated in Appendix Bl of iGravl (|l992l ). which is 
again valid over the temperature range 4000 < Tes < 7000 
K. 



Rq 



-3.925 X 10" 



'(Teff)'' + 8.3909 X 10"'° (Teff )^ 



-6.555 X 10"^(rcff f + 0.02245(reff ) - 27.9788. 



The stellar mass is then estimated via the main-sequence 
mass-radius relationship 



Mo 



'R^ 



1/0.8 



(B2) 



7000 - ^ 



a: eOOO 




ma 



Figure Bl. J — H colour index against temperature for FGK 
dwarfs from the Tycho-2 catalogue (only 1/15 datapoints are 
shown for clarity). 
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